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I. INTRODUCTION 

There is now a huge interest in the novel optical prop- 
erties of coherently-prepared atomic media, especially 
electromagnetically induced transparency (EIT) where a 
resonant coupling field (or pump field) coherently pre- 
pares an atomic sample to allow dissipation-free propaga- 
tion of a weak probe field accompanied by strong disper- 
sion. There have been applications of EIT to enhanced 
non-linear optical processes JjJ], laser cooling g], quan- 
tum non-demolition measurements and gain without 
inversion [Q . EIT is also the mechanism underlying the 
recent experiments in ultra-slow Q and ultra-fast group 
velocities || and is potentially of use for storage and re- 
trieval of quantum information using robust ground state 
coherences . There is also a continuing interest in ex- 
ploiting the strong optical non-linearities for the control 
of light by light at the single photon level || . A review 
of EIT and its applications can be found in ref |)| . 

The study of transient excitation of three level systems 
is a mature field, but there is still a need for experimental 
verifications of theoretical work. An early theoretical pa- 
per by Berman and Salomaa pp[ | compared the dressed- 
atom and bare-atom pictures, and considered transients 
after probe turn-on. Related dressed-atom transients 
for two level atoms are presented theoretically by Lu and 
Berman Jll]]. Theoretical three- level transient studies 
considering intial conditions have also appeared in Lu et 
al. [^2|. Harris and Luo |ll| studied transient EIT in 
the context of the energy required for the preparation 
of EIT. Li and Xiao |14| investigated the time required 
for the onset of EIT and Zhu considered the conditions 
required for observing inversionless gain in the transient 
regime for the V |l5| and A 1 16 schemes. 

Experimental work looking at dressed state transients 



includes phase shifting measurements in a two level sys- 
tem IE 71, fluorescence measurements on a three level A 



system jlq] and pump-probe experiments on a three level 
ladder system |19(| . Transient gain was first observed ex- 
perimentally in a sodium sample in a A configuration by 
Fry et al. pc| • This is also the only previous experiment 
we know of which analyzed transient dynamics associated 
with the turn-off of the coupling field. RF experiments 
on the N-V centre of diamond appear in plj . In our 
earlier work, we have demonstrated transient EIT J^] 
and Rabi oscillations and gain without inversion in 
a cold rubidium A system after rapidly switching on a 
resonant coupling field. 

The work presented here differs from all previous stud- 
ies in several ways. It extends previous work in tran- 
sient EIT by studying turn-on and turn-off transients, 
and the region of parameter space investigated is much 
broader than previously studied. This is the only work we 
know of which studies turn-off transients for off-resonant 
probe and coupling fields, inversionless gain away from 
resonance and transient dressed state interference away 
from either bare or dressed state resonance. Furthermore 
we present new analytical results to describe our exper- 
imental transients, derived using the Laplace transform 
method, so extending the two level approach presented 
in [^4| . Our theoretical results also generalize those of Li 
and Xiao Jl4j for mutual resonance to include two-photon 
dephasing and arbitrary initial ground state populations, 
and we predict a new frequency of oscillation in the probe 
transient response. We also present an analytical solu- 
tion, to first order in the probe strength, for the probe 
absorption when the coupling field is non-adiabatically 
turned off for arbitrary detunings of coupling and probe 
fields. Our analysis is performed explicitly in the bare 
state basis. Our results provide insight into how atomic 
systems become dressed by intense laser fields. 
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A natural classification of E1T response with switched 
fields is based on the relative time scales involved. In 
the adiabatic regime the switching is assumed to occur 
sufficiently slowly that the system evolves smoothly from 
one steady state to another. For EIT systems where 
the Rabi frequency is comparable to the atomic lifetime, 
the switching is adiabatic when it occurs on a time scale 
that is long compared with the relevant optical pump- 
ing times. For example, steady state A-type EIT with 
resonant fields is characterized by the ground state pop- 
ulations being in the non-absorbing superposition of the 
two ground states. If the intensity ratio or relative phase 
of the two optical fields changes, then the composition of 
the non- absorbing superposition state changes, and so 
the populations have to be pumped into the new non- 
absorbing state to maintain transparency. Thus the 
adiabatic condition requires the switching of the optical 
fields to be slow compared with the time taken for optical 
pumping from the absorbing state to the non-absorbing 
state. This is the regime used in the recent studies of 
optical information storage and retrieval Q], the single 
photon switch (2^] and its classical precursor [p6| . 

There is another adiabatic regime in pulsed systems 
where evolution occurs between ground states without 
spontaneous dissipation. In the work by Harris and 
Luo jl3| on the conditions for preparing EIT, this was 
achieved by using extremely large Rabi frequencies so 
that all system evolution to the dark state was achieved 
on timescales much shorter than the spontaneous emis- 
sion time. Not requiring such extreme Rabi frequen- 
cies, adiabatic switching using the well-known STIRAP 
schemes J2?j decouples the evolution from the sponta- 
neous emission, and so adiabatic evolution occurs on 
timescales independent of spontaneous emission time. 
The picosecond pulse experiment of Nottleman et al. |2f| 
is also effectively in the adiabatic regime. In their experi- 
ment the relative phase of the two Zeeman ground states 
oscillates by RF Lamor precession in a magnetic field, 
and the probe pulses are timed to arrive at the partic- 
ular times in the RF cycle when the ground states are 
in the non-absorbing state. We are not concerned with 
such regimes in this work. 

In this paper we are concerned with the non- adiabatic 
regime, where the switching time is very fast on the scale 
of optical pumping times, and all other relevant time 
scales, with the coupling field strength comparable to the 
spontaneous emission rate. The theoretical studies by 
Li and Xiao |l4| and Zhu jf(| assume (as we do) instan- 
taneous switching and are therefore in the non-adiabatic 
regime, as are the three experimental observations in A 
systems [^0[ ^3) where the coupling field was rapidly 
switched using a Pockels cell. These studies were re- 
stricted to special cases of resonant fields and some spe- 
cific detunings. 

There has also been interest in a further non-adiabatic 
time regime where effects are studied on timescales com- 
parable with the transition period, see for example ]29| ] 
where the Rabi frequency of the coupling field was com- 
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FIG. 1: (a) Model energy level diagram. The atomic energy 
levels are labelled \a), \b), and |c) in order of increasing energy. 
The \a) — c) transition is driven by a strong coupling field, C, 
with detuning, Ai and Rabi frequency, fii (t). The \b) — c) 
transition is monitored by a weak probe ,P, of detuning, A2 
and Rabi frequency Q.2 < T. (b) Experimental realization 
in 87 Rb MOT. The trapping fields T are detuned from the 
5S 1/2 F = 1 to 5P S/2 F = 3 transition by A T = —13 MHz and 
a repumper, R, is applied resonantly with the 5S1/2F = 1 to 
SPi/zF = 2 transition. The coupling field is applied to the 
5S1/2F = 1 to 5P 3 /2-F = 2 transition and is switched non- 



adiabatically using a Pockels cell. The probe is applied to the 
5S1/2F — 2 to 5P3/2F = 2 transition. Both coupling and 
probe fields were held at constant detunings for data collec- 
tion. 



parable to the transition frequency in a Galium- Arsenide 
semiconductor, and in the N-V centre of diamond |3(| 
where RF transitions at the Rabi frequency were further 
driven by intense fields. Such regimes are quite different 
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from those considered here as they go beyond conven- 
tional Bloch analysis of the atomic evolution. 

The organization of this paper is as follows. Sec- 
tion |n] outlines the Laplace transform method for calcu- 
lating transient evolution of an ideal A system, presents 
some simple analytical results for coupling field turn-on 
and turn-off, and gives an overview of the general predic- 
tions for arbitrary detunings of both fields. Section [II 
describes the experimental realization of a A system in 
laser-cooled rubidium and presents the experimental de- 
tails. The results of the turn-on and turn-off experiments 
are presented in Section [V together with the theoretical 
predictions. The results are discussed in Section [v|. 



II. TRANSIENTS IN A A SYSTEM USING 
LAPLACE TRANSFORMS 

We consider an ideal closed A system of three non- 
degenerate levels excited by monochromatic laser beams, 
as shown in FIG. 0(a). The sample is assumed to be 
optically thin and so propagation effects are not consid- 
ered. The atomic energy levels are labelled in order of 



increasing energy \a), \b) and |c). The \a) — |c) (\b) — |c)) 
transition is coupled by an optical field with frequency 
uj\ {0J2) and detuning Ai = uj\ — uj ca (A2 — UJ2 — cu c b)- 
The transition frequencies are defined u>p a — uip — Lu a , 
a, f3 — a, 6, c and in our system tOb a "C u ca , to c b- For con- 
venience we also define A21 = A2 — Ai. The state |c) is 
an excited state which decays by spontaneous emission 
to the ground (or metastable) states \a) and \b) at rates 
r cQ and T c b respectively. There is assumed to be decay of 
the (two-photon) coherence on the |a) — \b) transition at 
a rate r ba , but no incoherent coupling between these two 
states. The Rabi frequency of field 1 (2) is £!i = d ca • Ei 
(£^2 = d c fc • E2) where d a is the electric-dipole moment 
of the |a) — |/3) transition, Ej is the electric vector of op- 
tical field j, and we have chosen units such that h = 1 so 
that energies are measured in units of frequency. Opti- 
cal field 1 is the pump or coupling field (C) and field 2 is 
the probe (P) . For the purposes of this work we usually 
assume that the probe is weak, i.e. << ^cb although 
the density matrix equations are presented without this 
assumption. 

In the bare state basis, the density matrix equations 
of motion for the system in the rotating frame are 
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The method for switching between bare and dressed bases 
is given in JlOj ] - Because we measure probe absorption, 
which is proportional to Im[p bc (t)], we concentrate on 
this component of the density matrix. 

The standard approach to solving Eqs. (|l|) is to nu- 
merically integrate them (see for example |16|). This 
method has the advantage of conceptual ease and has 
been used to generate some of the theoretical turn-on re- 
sults presented in this work. However this approach 
has the deficiencies that it is not analytic and there- 
fore gives minimal insight into the underlying dynamics 
of the problem, and that it can suffer cumulative nu- 



(1) 

I 

merical errors. Accordingly, we have also considered a 
Laplace transform solution which can avoid these difficul- 
ties. The Laplace method (in appropriate limits) is more 
amenable to analysis and in this work we present what 
we believe to be several new analytical results, especially 
the turn-off equations. We find an extra frequency of 
oscillation when the system responds to turning on the 
coupling field, and suggest the regime where this might 
be observable. Even when the Laplace method cannot 
easily produce simple explicit formulae, Eqs. (0) are lin- 
ear so that any solution depends only linearly on the ini- 
tial conditions, and so avoids the problem of cumulative 
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numerical errors, making the solutions more robust. 
Any solution to Eqs. (|l|) has the form 

p a p{t)=Y J W hlt - (2) 
i 

Our goal therefore is to determine coefficients ai and 
The difficulty lies in the fact that the coefficients depend 
on some or all of the system parameters, including f2i 
and fl 2 . 

Defining the Laplace transform of p a p (t) to be r a p (p) 



we have H| 




Pa/3 (t) 



and the transform of p a p {t) is pr a /3 (p) — p a a where 
Pap — P a P (0)- Then Eqs. (|]) generate nine coupled 
algebraic equations for the r a p (p) . By manipulating 
these equations, one finds four closed equations for r ac , 
T"ca, fbc, and r cb . One of them is 



(Q 2 /2y 



P - iA 21 + r 



bo 



+ r» W [» (P) " 1] - ric (p) 12i5i s (p) + r c , (p) 1 f!l!i2 



p 4 



p 4 p 4 

o r 



P° cc s(p)} 



p - iA 2 i + T ba 

n 2 



Pab 



2 p - iA 2 i + r ba 



(3) 



r 



with 



where 



s(p) 



P 



P + {Tea + r 



cb ) 



The remaining three equations can be obtained from 
Eq. (||) by (i) taking the complex conjugate of Eq. (||), 
treating p as real, (ii) interchanging all labels a and fo, 
and 1 and 2, in Eq. (0) and (iii) taking the complex con- 
jugate of the latter 'interchanged' equation. Keeping p 
real in these procedures is merely a formal device. Once 
the equations are written by the scheme given, p may 
take on complex values, as it must when one effects the 
Laplace inversion to get p a p (t) . 

The four coupled equations allow for arbitrary initial 
values p° o. Finally when r ac , r bc , r ca and r cb are known, 
one has 



"fab (j>) 

n a ip) 



Pab + i^at ( P ) 



ac (p) 



p - iA 21 + r 6a 

i^r ca (p) -i^r bc (p) 



Pba + * 



p + iA 21 + T b 



(4) 
(5) 



and the diagonal elements are 



1 



Taa (p) = - [p°aa + ^caTcc (p)] 
p 

+*--7T l r ca (p) - r ac (p)} 
P 2 

Ub ip) = - [p bb + T cb r cc (p)] 



-*--7T [ r cfc (?) ~ r bc ip)] , 

P 2 



(6) 



(7) 



■ip) 



ip) 



ip)\ 



[rbc ip) 



+ r cb 

Tab ip)] 



(8) 



Having solved the equations for the r a p (p), one generates 
the p a p (t) by Laplace inversion, namely 



Pa/3 {t) 



2ni 



dpe pt r a p (p) , 



(9) 



where the contour C in the complex p-plane runs verti- 
cally from a — zoo to a + ioo and a is real and chosen 
sufficiently large so that C runs to the right of all poles 
of r af3 (p). 

Equations (H) to (H) yield solutions of the form 



r a ip) = 



Cap (P) 
VapW 



where the C a p (p) are polynomials which depend on the 
parameters of the problem, and in particular, linearly 
upon one or more of the initial values /?% (0). The 
Tap ip) are polynomials too, but do not depend on the 
initial values. Generally the zeros of V a p (p) in the com- 
plex plane are distinct, so that r a p (p) may be expanded 
as partial fractions 



r aP ip) = 



ip) 



a/3,1 



(10) 
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where P a p.i is the I root of the equation V a p (p) = 
and the !Z a p.i (p) are determined using standard algebra. 
Using this in the contour integral, Eq. (||) gives 



Pap 



(11) 



Although the principles of this approach are simple, the 
difficulty lies in finding the zeros of the V a p (p) , which 
may have to be done numerically or by some approxima- 
tion. Most roots P a p,i give rise to damped oscillatory 
terms or to simple damping, but should there be a pole 
at the origin, its contribution does not decay and repre- 
sents a stationary state. Such a long-time limit, should 
it exist, is 



lim p a p (t) = TZ a p (0) = lim 

t— >oo p^O 



pCg/3 (g) 

V a p (p) 



(12) 



To facilitate matters, at very little cost, we shall hence- 
forth take r ca = r c f, = r. 

We now present a simple illustrative example of the 
Laplace transform method, showing explicit expressions 
for optical pumping due to the probe only. This is fol- 
lowed by application of the method to the transients fol- 
lowing turn-on and turn-off of the coupling field. The 



J 



latter solutions are provided to first order in the probe 
intensity. For simplicity the turn-on expression is only 
given for resonant fields, however the turn-off result is 
presented for arbitrary detunings. This section ends 
with a graphical overview and discussion of the general 
turn-on and turn-off dynamics. 



A. Optical pumping by a probe 



As an illustration, we consider optical pumping by 
the probe in the absence of a coupling field (Sli = 0). 
This simple illustrative example encapsulates the Laplace 
method. A Schrodinger equation approach to this famil- 
iar problem appears in |32j] . We suppose that when the 
probe is turned on, at time zero, all three states may be 
populated but that there are no coherences, i.e. p Q a ^ may 
be non-zero for a = (3 only. Then, when Eq. (|^) and the 
associated interchanged equations are solved, one finds 
that only p bci p cb , p aai p bb and p cc (= 1 - p aa are 
in general non-zero. Use of Mathematica 4.1 |33| to do 
the algebra, yields 



n c (p) = 

r aa (p) = 

nbip) = 



^2 j£ (1 - p° aa + P ° bb ) + P (2pg b + PL - 1)] [i (P + r) + Ag] 
2 Dip) 

p [r + Paa (r +p) - P ° bb r] [Ai + (p + rfl + 2 ( P + r) (r + 2p° Q p) (n 2 /2) 



P D(p) 

[r (1 - pL + p° bb ) + P° bb p] [a| + (p + r) 2 l + 2 (p + r) ( P ° aa - 1) (n 2 /2) 



D(p) 



(13) 
(14) 
(15) 



where 



£(p)=p(p + 2r) A 2 2 + (p + r) 2 +2( P + F)(2p + r)(n 2 /2) 



Furthermore, r cb (p) is obtained by complex conjugation 
of r bc , holding p real, and then allowing p to be complex. 
Then the p a p (t) follow by Eq. @. Applying the limit in 
Eq. ([l2]) to these quantities shows that all p a p (t) except 
p aa approach zero as t — > oo, and p aa {t) — > 1 in that 
limit. This simply means that even a weak probe field 
will pump the system into state \a) after a long time, 
the usual result. A close look, however, shows that for a 
weak probe field (0 2 < r) there are three important time 
periods. Since D (p) is a fourth degree polynomial in p, 
it has four roots. The real parts of three of them are 
of order — T or — 2r but the fourth root is much smaller. 
To see this, we set D (p) = and assume that terms of 



order p 2 and higher are small. The fourth root P4 is then 
approximated by 



Pi 



A 2 +T 2 



Thus, there is an initial period of damped oscillations in 
Pap (t), which die out in times of order or (2r) 
Optical coherences form in an intermediate period after 
the initial dampening, and then a longer period obtains 
when p bc it), p bb (t) and p cc (t) very slowly damp to zero 
and p aa it) slowly approaches unity. This latter time 
period derives from the probe optical pumping. For 
small f2 2 , the intermediate period is long and lies between 
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1/ (2r) and l/p 4 - For this epoch, Eqs. (Ill),® and (£|) 
are well approximated by expanding them to first order 
in CI2 ■ When this is done and the limit |l2| is taken, one 
finds 



Pbc (t) 
Paa (t) 
Pbb (t) 



4 A 2 - iT ' 

1 + Paa ~ Pbb 

2 

1 - Paa + Pbb 



More generally, although the equation D (p) = is of 
fourth order and exactly solvable, numerical values for 
Pa/3 (t) are perhaps best achieved by numerically finding 
the roots p\ to p^. 



B. Turn-on transient with resonant fields 



We now shift the time origin to the instant of turn-on of 
the coupling field and use Eqs. (|J) to (|J) to find Im [pb c (i)] 
after t = 0. It is assumed that the probe has been on 
for some time in order to establish values for p® a , p bb and 
Pbc = (Pcb) * ■ F rom t ne previous discussion we know that 
in our case, these will be the only non-zero values of p a g. 
The probe is assumed weak, VL 2 <C T; therefore we treat 
the problem to first order in the probe intensity. This 
corresponds to the experimental situation to be described 
in Sections III and IV. For simplicity we here assume 
Tb a = and resonant excitation, i.e. Ai = A2 = 0. We 
then find, |35|] , that the Laplace transform of Im [p bc (t)] 
(to first order in O2 but all orders in Oi ) follows from 



J 



r bc ip) - r cb (p) 
2i 



Im I 



( P ° aa - P l b ) n 2 /2 



n 2 



rp° aa + (p bb -i)(p-r) + 2p° aa p 



P T+p 2 + (fii/2) 2 pT + p 2 + (fti/2r 2 P (P + T) (p + 2T) + 2 (2p + T) (fti/2) 



(16) 



r 



It is worth looking at the partial contribution given by 
the first term on the right hand side of Eq. (TLa) , namely 



im(p° bc ) P / P r + p 2 + (n 1 /2y 



for which the poles are 



p = --±(i/2)^q 2 -t 2 . 

The inverse transform is easily found to give 



lm[p bc (t)} = Im^Je-™/ 2 . 



ft 



-l m (pl)e-^sin^ 



(17) 



where / = — T 2 . A similar result is given by Li 

and Xiao jl4| and corresponds to standard damped Rabi 
nutation. 

More generally, the contributions of p {) aa and p bb may 
not be ignorable, especially if Im (p bc ) is small or zero. 
Therefore we look at the full expression, Eq. ([l6|), for 
which the first two terms are easily inverted. The de- 
nominator of the third term, however, is a cubic poly- 
nomial in p and can be somewhat complicated. But as 
we are primarily interested in the case when the coupling 



field is strong, namely when £!i ^S> T, then the roots, p\ 
and pj. in the first two terms on the right hand side of 
Eq. (Jlq) can be estimated as 



T Mi 
and those of the third term as 



P3 
Pi 



r 

"2 5 



5r 



neglecting terms of order T /f2i. Note that pi and p 2 give 
rise to damped oscillations of angular frequency fii/2 as 
seen from Eq. (^), that p^ yields a simple damped term, 
and that P4 and p§ give rise to damped oscillations of 
angular frequency f2i. Using these roots, it is a mat- 
ter of algebra to derive the following approximation to 
Im [p bc (t)}: 



Im [p bc (t)} = $12 (t) + $3 (*) + * 45 (*) ; 



(18) 



where 
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$12 (*) 
$3 (t) 

$45 (t) 



-rt/2 



ft-, 12(1 



■A) 



r 9 + (40 x /rr 

^2 e _ 5rt / 4 r 12 fi 



nit 



P° bb ) 



9 + (40i/r)" 



r 



cos (filt) 



r 



P ° bb , 



P ° bb 



r 



~2~ 



2pL r 



Oi 



sin (Hit) 



(19) 
(20) 

(21) 



We have compared p3j the exact numerical inversion of 
Eq. (Rq) with the approximations in Eqs. (19) to (21). 



The agreement is qualitatively good even for values of 
f2i/r greater than 2. 

In Eq. (fl8|) ; $12 (t) is the contribution of the poles 
at pi and p 2 - That part of $12 (i) which is propor- 
tional to Im (Ph c ) reduces in the limit of f2i 3> T to 
Eq. (17]). $3 (t) results from the pole at p% and is 
strictly monotonically decreasing with time. $45 (t) re- 
sults from poles p^ and p$ and oscillates (in the limit 
ill S> r) twice as fast as $12 (i). In all the transient 
experimental and theoretical studies we know (including 
this present experimental study), the system behavior is 
dominated by $12 (t); however, should Im (p° c ) be small 
or zero with approximately equal ground state popula- 
tions (i.e. p° aa ~ p® b ), then the faster oscillations of <f> 45 
might be observable. One way to do this would be to 
switch on both the probe and pump non-adiabatically at 



t = 0. In this case we would have (in general) non-zero 
Paa, Pbb and Pec = 1 -Paa-P bb , but vanishing off-diagonal 
components. A discussion of the faster oscillations using 
the 3D Vector Model Jjl |j§ appears in the Appendix. 

Equations @ to Q were derived putting r;, a = 0. 
In experimental realizations, a non-zero value of Tf, a can 
result from collisions, trap inhomogenieties and other ef- 
fects. We can generalize Eq. ([l7]) including this dephas- 
ing, to yield 



Im [p bc (t)} 



Im( P ye-*( r + r -)/ 2 ^-^sin( 



ft 
2 



+ Im(p° c ) e -'( r + r -)/ 2 cos( if), (22) 



where /' = y Q\ — (T — Tba) ■ We can also obtain the 
general forms of Eqs. (H) to @ (for Q x > T): 



J 



$12 (t) 

$3 (t) 
$45 (t) 



- p -t(r+r,, a )/2 



Im(p° c ) 



-t(r+r 6o )/2 

ni 
n 2 



^2 -SCt/ 4, 



« 2 r bo 5p° 6 + p° a -2 



Im (p 6c 



r-r 

n 5 



of 

raa 



COS 



r fcQ - I (1 - P ° bb ) (r - 4r bQ ) e - rt / 2 
3n~p° b r ba i-2p| 



aa 



"2" 



1 - 2n° - n° 
cos (firf) + ^ — ^ sin 



r 



C. Turn-off transient for arbitrary detunings of 
pump and probe 

The turn-off case is relatively simple as only the probe 
field acts during the transient period. Thus we are con- 
cerned with optical pumping by the probe, as considered 
earlier, but in this case the initial values arc different 



and we shall work to first order in the probe field ■ 
We suppose that both the coupling and probe fields have 
been turned on for a time long compared to T _1 and the 
coupling field is then switched off at a time now taken to 
be t — 0. We find |33j that the inverse Laplace transform 
of p bc (t) is 



J 



r bc (p) 



Pi 



+ i 



n 2 (p - r) (1 



p° ) 

raa } 



pg 6 (2p + r) 



p(p + 2T) 



o 



r 
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The initial values p° aa: p bb and p° are found by solving for the r a p (p) during the epoch when both fields act, and then 
taking the long-time limits, Eq. (|12|). Thus we have 



Paa 



Pi 



Pi 



o 



f2; 



1 + 



A 2 i - iT ba 



- a 2 a 21 + r ba r) + i (A 21 r + A 2 r ba ) 



so that, 



r&c (p) = 



p + r + iA 2 

This expression contains simple poles at p = and p = — T — iA 2 so that for t > one has 



„— ri„— iA 2 t 



2 \ 1 + zA 2 

Finally, using Eq. ( p3|) for p° c , and taking the imaginary part of p bc (t), gives to order f2 2 /r 



Im [pfe (t)] = - 



-rt 



[A 2 sin (A 2 t) -Tcos (A 2 <)] 



2 

2 



r 2 + A 2 



-r< 



[ta 2 21 + r ba (rr ba + nj/4)] cos (A 2 t) 
(ir ba - a 2 a 2 i + «2/4) 2 + ( A 21 r + A 2 r ba ) 2 



^2 rt 
2 



A21 



%-A 2 (A 21 +r 2 a /r) sin(A 2 i) 



(rr ba - a 2 a 21 + nyiy + (A 21 r + A 2 r ba r 
1 



(23) 



(24) 



(25) 



which is the required turn-off transient. It is interesting 
to note that r ba enters Eqs. (|24|) and (p^) only through 
the initial conditions. This may be understood by re- 
calling that r ba refers to a two photon dephasing, and so 



will not be dynamically important when only the probe 
field is acting. We note that, at t = (the instant of 
turn-off) one has, from Eqs. (|23|) or (pj 



Im(pg c ) = 



n 2 



TA 2 21 + r ba (rr ba + nj/4) 



2 (rr ba - a 2 a 21 + nl/4) 2 + (A 21 r + A 2 r ba ) 2 



which is the general expression for the steady-state EIT 
lincshape for arbitrary detunings. 



D. Overview of transient response 

We now present an overview of the turn-on and turn- 
off transients displayed as pseudocolour plots showing 
Im [p bc (t)] as a function of time t and probe detuning A 2 . 
The parameters are chosen to correspond closely to the 



experimental conditions described in detail in the follow- 
ing sections. We use the spontaneous radiative decay rate 
(2r = 5.68 MHz) and the effects of uncoupled absorptions 
are not included in the model. (Uncoupled absorptions 
are extra absorptions which derive from degenerate en- 
ergy levels: see for example (36) for a full description). 
In order to generate the pseudocolour plots, we numeri- 
cally turned on the probe field at a finite negative time 
(t = — 10r) which was sufficient to build up a steady state 
Pfc. Initial populations (at t = — 10T) were chosen so 
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FIG. 2: Pseudocolour plots of Im [pb c (t)]as a function of time, t, and probe detuning, A2 for (a) turn-on; resonant coupling 
field: Ai = MHz, (b) turn-on; detuned coupling field: Ai = —23 MHz (c) turn-off; resonant coupling field: Ai = MHz and 
(d) turn-off; detuned coupling field: Ai = -23 MHz. In all cases we have used: fii = 45 MHz, tt 2 = 1 MHz, V = 5.68 MHz and 
Tba = 3.4 MHz. (The effect of uncoupled absorptions are not included in these plots.) 



that the populations at time t = corresponded closely 
to the experimental conditions (p aa (t = — 10T) = 0.05, 
Pbb (t = — 10F) = 0.95). We term this the 'linear regime' 
as small variations will only introduce extra global scal- 
ing factors to the amplitude of the probe response. 

The turn-on transients are shown in FIGS, ^(a) and 
(b), where the coupling field is switched on at t = 0. 



We consider first the case of a resonant coupling field, 
FIG. ||(a). Scanning the entire figure from left to right, 
we first see (for t < 0) the Lorentzian probe absorption 
feature of spectral width 2r as the probe drives the bare 
\b) — |c) transition. Then in the transient regime from 
t = to several we find modified Rabi nutations 

due primarily to the coupling of the ground state \b) to 
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the dressed states |±) = (\a) ± |c)) /\/2- If the probe 
were to couple |6) to a single dressed state, say |+), it 
would undergo decaying Rabi oscillations at the detuned 

frequency J (A2 + Oi/2) 2 + f2|. The peaks of these os- 
cillations would lie on the curves A 2 = ±2nn/t + H.i/2 
for t <C 2n/Q,2 and n — 1,2, ... . Oscillations of this 
kind are clearly evident in FIG. |(a) for each of the 
dressed states. However, the fact that the probe couples 
I b) to both dressed states simultaneously evidently, re- 
sults in interference between the Rabi oscillations due to 
the probe interacting with each dressed state separately. 
This suggests that the central pattern in FIG. |^(a) repre- 
sents transient dressed-state interferences. In addition, 
for A 2 « 0, the preparation of the system before the cou- 
pling field is switched on induces relatively strong absorp- 
tion and a larger value of Im(pf, c ) at t = 0. This leads to 
enhanced Rabi nutations along the A 2 = line for t > 0. 
Finally, for long times, we see the two well-resolved ab- 
sorption features, each of width T, corresponding to the 
well known Autler-Townes doublet |3lJ. We note that 
our pseudocolour plots show a superficial similarity to 
the three dimensional plots presented by Lu et al JlJ] 
but the different initial conditions between the two cases 
studied result in quite different dynamics. 

We now focus on the case where the probe too is fixed 
on resonance, i.e. the time-line A 2 = in FIG. ||(a). 
Here, for t > 0, we see gai n and absorption cycles of the 
kind discussed in refs ju], |l6| and observed by us pre- 
viously |^2|, ^3). For a strong coupling field (f^i/r ^> 1) 
this transient (A 2 = 0) is well described by Eq. (22). 
Two other lines of note are at A 2 = ±J7i/2 where trans- 
parency for t < changes transiently to absorption as 
the probe begins to resonantly monitor a transition to a 
dressed state. 

Figure |2|(b) shows the case of a detuned coupling field. 
Again, we observe a transient regime featuring absorp- 
tion and gain regions evolving to the resolved dressed 
state absorption features centered at A 2 = Ai/2 ± 
v/(fif + A 2 )/4, but in this case we can clearly distin- 
guish between the dominant or major dressed state and 
the narrower minor dressed state. In the bare state pic- 
ture, in the limit of large Ai these two features are associ- 
ated, respectively, with one-photon absorption (\b) — |c)) 
and two-photon absorption (|&) — |c) — \a)). 

We now consider the turn-off transients shown in 
FIGS, ^(c) and (d), where both fields are assumed to 
have been on for all negative times to establish the ini- 
tial conditions, and the coupling field is switched off at 
t = 0. The dominant feature in the transient regime is 
the 1/t broadened absorption region where the probe be- 
gins to drive the bare transition. Of more interest are 
the strong transient gain pea ks that can be seen along 
the two lines A 2 = Ai/2 ± + A 2 )/4) at the levels 

associated with the two dressed state transitions before 
turn-off. This is simply the off-resonant probe moni- 
toring the relaxation of a two-level atom with specially 
prepared initial conditions. We note that these tran- 
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FIG. 3: Schematic of experimental beam configuration. The 
coupling, C, and probe, P, beams make an angle of about 20° . 
PD is an avalanche photodiode, PC is a Pockels Cell and Pol 
a polariser. P is linearly polarised in the horizontal plane 
while C is linearly polarised in the vertical plane when the 
PC is on and linearly polarised in the horizontal plane when 
the PC is off. 



sient gain regions are tunable by varying the detuning 
and strength of the coupling beam, and are remote from 
any uncoupled absorptions in a real A system that might 
exist in the A 2 = region |36[ Q . We note also that the 
entire (t, A 2 ) region of this figure is described by Eq. J25|). 



III. EXPERIMENTAL SET-UP AND 
PROCEDURE 

The A system was realized in a Rb 87 MOT substan- 
tially similar to the one described in [ p3| . The sample 
contained about 5 x 10 7 atoms in a region of diameter 
about 3 mm. The relevant energy levels and excitation 
scheme are shown in FIG. 0(b). A schematic show- 
ing the layout of coupling and probe beams is shown 
in FIG. The coupling beam had an average in- 

tensity 100mW/cm 2 with a roughly elliptical pro- 
file 2 mm x 4 mm. The probe beam had an average in- 
tensity w 0.03 mW / cm 2 with circular profile of diam- 
eter about 1 mm. The beams propagated in the sam- 
ple with orthogonal linear polarizations and with an an- 
gle of about 20° between their directions of propagation. 
All laser beams were derived from external-cavity laser 
diodes, using master-slave injection where appropriate. 
Non-adiabatic switching of the coupling field was realized 
using a Linos Photonics LM10 Pockels cell with a cus- 
tom designed high voltage switch. The rise time of the 
switch and associated detection circuitry was less than 
6 ns. The excited 5P 3 / 2 state lifetime is 28 ns and the 
relevant Rabi periods are of similar magnitude and so the 
non-adiabatic condition is fairly well satisfied. Because 
the fall time of the switch is very long, the Pockels cell 
could only be operated in either fast turn-on or fast turn- 
off mode, necessitating separate experiments to observe 
turn-on and turn-off transients. Changing between turn- 
on and turn-off experiments was realized by rotating the 
plane of polarization of the coupling field by 90° before 
entering the Pockels cell. 

The probe detuning was set by varying the cavity 
size of an external-cavity grating-controlled diode laser 
(ECDL) using the voltage offset of the piezo-electric 
transducer, which was monitored on a voltmeter. This 
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voltage was calibrated using the saturated absorption 
spectrum in a standard Rb cell. Drifts in the cavity 
size due to temperature fluctuations were checked and 
corrected for by determining the voltage which put the 
probe field on resonance at the beginning of each tran- 
sient experiment. Nevertheless, the absolute detuning of 
the probe field could not be determined to an accuracy 
better than about half a linewidth (as 3 MHz), except 
when the field was locked to resonance. Standard detun- 
ings used were generally much larger than 3 MHz, with 
dressed state detunings being as large as 36 MHz. The 
measured transients were stored in a Tektronix TDS520B 
digital oscilloscope. Instabilities in our control method 
limited the maximum number of averages to about 50. 
Because of small drifts in the intensities of our laser 
diodes (which were run in constant current mode for max- 
imum frequency and mode stability) , probe transmission 
levels were measured with the MOT turned off between 
transient runs and the transmission level at the time of 
the experiment inferred using a linear regression. 

Our experiments were carried out in the optically thin 
regime with maximum probe absorption of about 15% at 
resonance. The probe Rabi frequency of about 4 MHz 
gave a good signal to noise ratio without producing any 
measurable power broadening. Under these conditions 
the probe absorption increases linearly with probe inten- 
sity and it was found convenient to normalize our probe 
absorption/transmission levels to the steady state ab- 
sorption in the absence of the coupling field. We use 
the same method as that described in []23j where we de- 
fine the maximum probe absorption in the absence of 
a coupling field as the zero transmission level. Thus, 
if we define the voltage measured on our photodetector 
as V, the signal for maximum probe absorption (mini- 
mum probe transmission) V m i n and the signal without 
any probe absorption Vo, then our scaled transmission 
levels for a detected signal are defined by 



Vo — Vmin 

This definition for the transmission has the advantage of 
being simple and robust against day-to-day fluctuations 
in the number of atoms trapped in our sample (provided 
we are in a linear absorption regime). It is also easy to 
make direct comparisons with theory by noting that any 
detected changes in signal level will be proportional to 
Im(pt,c), with values of T > 1 associated with gain and 
T < 1 indicating absorption. 

With on-resonant or near-resonant probe fields our 
experimental transients are accompanied by the time- 
independent effects of uncoupled absorptions. These 
have been fully described in our previous paper |3(| . For 
the purposes of this work it suffices to treat the system 
as an ensemble of identical three level atoms in an ideal 
A configuration, as in Section I, together with a similar 
ensemble that does not interact with the coupling field. 
In our turn-on experiments, uncoupled absorptions were 
found to contribute about 20% of the maximum reso- 



nant absorption and are assumed to have a Lorentzian 
line shape centered at zero probe detuning. 



IV. EXPERIMENTAL RESULTS AND 
DISCUSSIONS 

We first present experimental steady-state probe fre- 
quency spectra (EIT traces) for a resonant and an off- 
resonant coupling field with the same parameters that 
apply in the turn-on transient experiments. The central 
region of FIG. |^(a) shows the probe spectrum for the res- 
onant coupling field. The dressed state absorption peaks 
are clearly resolved on either side of the uncoupled ab- 
sorption at A2 = 0. The Rabi frequency of the coupling 
field can be inferred from this figure to be fli =45 MHz . 
This value is consistent with that obtained from beam 
parameters and absorption line strengths. Also seen are 
the absorption peaks on the 5Si/ 2 F — 2 to 5P 3 / 2 F = 1 
transition at A 2 = -157 MHz and the hSi/ 2 F = 2 
to 5P 3/2 F = 3 transition at A 2 = 267 MHz. These 
peaks were used for frequency calibration purposes. Fig- 
ure |^(b) shows the probe spectrum when the coupling 
field is detuned by Ai = —23 MHz. The major dressed 
state absorption peak is seen at A 2 = 14 MHz and the 
minor one at A 2 = —37 MHz. Uncoupled absorptions 
are barely resolved in the wing of the major dressed state. 
We note that the linewidths in these spectra are approx- 
imately 1 1 MHz, or nearly twice the spontaneous decay 
rate of the 5P 3 / 2 levels. We attribute this additional 
broadening to inhomogencitics in the sample principally 
due to light shifts caused by the coupling beam and the 
trapping beams with smaller contributions from laser 
linewidths. Accordingly we have taken a phenomeno- 
logical value for the half linewidth of T = 5.5 MHz in our 
theoretical transients rather than the standard value of 
2.8 MHz, as discussed in |3(| . The two-photon dephasing 
rate is taken to be r 6a = 0.6r ||. 

We now present experimental turn-on and turn-off 
transients for the resonant and the non-resonant coupling 
field, together with the corresponding theoretical curves. 
Turn-on theory transients were obtained by numerically 
solving Eqs. (jl|). Turn-off theory transients were ob- 
tained by fitting the experimental data to Eq. (|25|). In 
each case representative selections of probe detunings are 
chosen so that horizontal lines across FIGS. ||(a)-(d) are 
explored. 

When we compare the experimental and theoretical 
curves in FIGS. |||| we find excellent agreement in the 
transient behavior. We therefore use our theoretical re- 
sults to determine A 2 which we find to be consistent with 
detunings derived from the piezo-offset. It has been less 
easy to make direct comparisons of the overall transmis- 
sion levels and there are minor inconsistencies when com- 
paring experimental with theoretical levels in our tran- 
sients. We believe that these inconsistencies are due to 
difficulties in correcting for all of the intensity fluctua- 
tions and drifts in our lasers on the timescales of our 
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FIG. 4: Steady state probe transmission spectra, (a) On- 
resonance coupling field with Rabi frequency Q,\ = 45 MHz 
showing dressed state absorptions at detunings A 2 = ±f2i/2 
and uncoupled absorption at A2 = 0. (The peak on the 
5Si/2,F = 2 to 5Ps/2,F = 3 transition has been trun- 
cated to fit on the page). (b) Off-resonant coupling ex- 
citation with Rabi frequency fii = 45 MHz and detuning 
Ai = -23 MHz. The major dressed state is at A 2 = 14 MHz 
and the minor one at A2 = —37 MHz. Uncoupled absorp- 
tions are barely resolved in the wing of the major dressed 
state. In each case note the presence of extra absorption 
peaks corresponding to the 5Si/2, F — 2 to 5P3/2, F = 1 and 
5S 1/2 ,F = 2 to 5P 3/2 ,-F = 3 transitions at A 2 = -157MHz 
and A 2 = 267 MHz respectively for frequency calibration pur- 
poses. Each trace is an average over 50 scans. Day to day 
fluctuations in signal levels mean that the absorption scales 
in (a) and (b) are different. 



A. Turning on the coupling field 

For these experiments the coupling field is turned 
on at time zero after the probe field and MOT fields 
have been on for a long time. The left-hand traces in 
FIGS. ||(a)-(c) show the experimental results for turning 
on a resonant coupling field (Ai = 0) of Rabi frequency 
fii = 45 MHz. The right hand traces show the cor- 
responding theory curves, with the steady transmission 
levels adjusted to take account of the uncoupled absorp- 
tion. Figure |5|(a) shows the well-known case of zero probe 
detuning (A2 = 0) exhibiting significant gain (T = 1.35) 
despite the presence of uncoupled absorption. This curve 
is substantially similar to FIG. 3(b) of |§f]. The long- 
time transmission settles down to the value of T = 0.8 
showing the effects of uncoupled absorption as well as ef- 
fects due to the phenomenological two-photon dephasing 
rate. In FIG. ||(c) the probe detuning is 20 MHz and 
so the probe transmission is almost unity before turn-on. 
After turn-on, the transmission falls transiently as the 
probe begins to resonantly couple the dressed state. The 
long-time transmission is close to the value 0.5 which is 
expected from the equal absorption strengths of the two 
dressed states. We term this phenomenon, dressed state 
filling, because the increase in probe absorption and the 
associated movement of coherences and populations are 
accompanying the transient dressing of the \a) — |c) tran- 
sition. This has previously been called three-level free- 
induction decay |l(| . Figure ||(b) shows an intermediate 
case with probe detuning of A2 = 8 MHz. No significant 
differences were found between negative and positive de- 
tunings of the probe, as expected from the theory, and 
so negative detuning traces are not shown. 

We now consider the case where the coupling field is 
detuned from resonance by Ai = —23 MHz. Other pa- 
rameters and the dynamics before turn-on are the same 
as before. Experimental and theoretical transient results 
are presented in FIGS. |](a) to (e). The first three figures 
show (a) the case of a resonant probe, (b) the filling of 
the major dressed state (A2 = 15 MHz), and (c) the fill- 
ing of the minor dressed state (A2 = —37 MHz). Trace 
(d) (A2 = —6 MHz) shows the case where the probe is 
detuned just to the red of the bare state absorption line. 
It is interesting because of the clarity of the signal ob- 
tained coupled with the relatively high gain (T w 1.26) 
making this a promising region for further transient gain 
experiments. Finally in FIG. |(e) (A 2 = 9 MHz) the 
probe is detuned by a similar amount to the blue and 
shows enhanced absorption dips. 



B. Turning off the coupling field 



experiments, especially as they relate to mode hops in 
the probe laser as it was detuned. 



The experimental conditions for turn-off were different 
in some respects from those for turn-on. The maximum 
coupling-field intensity passing the Pockels cell for turn- 
off was fii = 46 MHz, slightly greater than that for turn- 
on; also the uncoupled absorption levels are slightly dif- 
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FIG. 5: Turn-on transients for varying probe detunings with a resonant coupling field: (a) A2 = OMHz (resonant probe), (b) 
A2 = 8 MHz, and (c) A2 = 20 MHz (dressed state filling). In all cases Qi — 45 MHz and 50 averages were taken. 



ferent, probably as a result of optical pumping amongst 
the Zeeman substates by the coupling field before turn- 
off. We have found that some of the experimental turn-off 
transients are remarkable for the high gain and clarity of 
signal obtained. We attribute the clearer signals to the 
fact that the system is simpler without the inhomoge- 
neous light shift effects of the coupling field. 

Using the analytical result, Eq. ( p5]) with the addition 
of the effects of uncoupled absorption, we have been able 
to perform fits to the experimental data using Origin 
6.0284 |^. The results of these fits are presented as 
the theoretical traces. The interdependence of many of 
the parameters mean that some parameters cannot be 
adequately determined using this method. Fortunately 
the probe detuning, A2, can be well determined using 



this method, and the agreement between experiment and 
theory shows that our system is well described by our 
three level analysis. Except where mentioned, the values 
for A 2 quoted below were determined using the curve 
fitting. 

The experimental and corresponding theoretical tran- 
sients for a resonant coupling field are presented in 
FIGS. 0(a) to (d). Here t = defines the turn-off time. 
Case (a) for a resonant probe shows the transmission 
dropping from the uncoupled absorption level of T rs 0.6 
to T — as the probe fills the bare state transition. 
Of more interest is case (b) where the probe is detuned 
to a level corresponding to a pre turn-off dressed state 
(A 2 = —22 MHz ~ -fii/2). We term this dressed state 
emptying by comparison with the dressed state filling re- 
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FIG. 6: Turn-on transients for varying probe detunings with an off-resonant coupling field: (a) A2 = OMHz (resonant probe), 
(b) A 2 = 15 MHz (major dressed state filling), (c) A2 = —36 MHz (minor dressed state filling) (d) A2 = —6 MHz and (e) 
A2 = 9 MHz. In all cases fii = 45 MHz, Ai = —23 MHz and 50 averages were taken. 
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FIG. 7: Turn-oil transients for a strong, on-resonant coupling field: (a) A2 = OMHz (bare state filling), (b) A2 = 22 MHz 
(dressed state emptying), (c) A2 = 17 MHz and (d) A2 = —36 MHz. In all cases Qi = 46 MHz before turn-off and the 
uncoupled absorption was found to be about 20%. 
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FIG. 8: Turn-off transients for a strong, off-resonant coupling field: (a) A2 = 0MHz (bare state filling), (b) A2 = 15 MHz 
(major dressed state emptying), (c) A2 = —40 MHz (minor dressed state emptying) and (d) A2 = —29 MHz. In all cases 
Ai = —23 MHz, fli = 46 MHz before turn-off and the uncoupled absorption was found to be about 20%. 
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FIG. 9: Close-up of experimental FIG. |^(c) showing the mi- 
nor dressed state emptying (solid curve) with added exponen- 
tial decays overlayed (dashed curves). The inverse lifetime 
of the decay was found to be 5.5 MHz which corresponds well 
with the measured probe spectral linewidths. 



ported above. The very clear signal features a gain peak 
rising to T = 1.35, making this another promising region 
for further study of transient gain. We note that this 
is a tunable gain region and is remote from any uncou- 
pled absorptions. The figure also shows two full periods 
of oscillation. The measured period of oscillation was 
found to be 46 ns which is consistent with the probe de- 
tuning of —22 MHz determined from the cavity off-set 
voltage. Figures 0(c) and (d) were taken for detun- 
ings of A2 = 17 MHz and A2 = —36 MHz respectively. 
These results show that cycles of gain and absorption can 
be present even when there is little or no change in the 
steady state transmission. 

The turn-off results for a detuned coupling field are 
shown in FIG. ||. Case (a) shows the transient filling 
of the bare state and is similar to FIG. 0(a). Cases 
(b) and (c) show the major and minor dressed states 
emptying respectively. As in the case of a resonant cou- 
pling field, these transients are characterized by ringing 
and gain. In case (b) (A2 = 15MHz, corresponding 
to the major dressed state) these effects are somewhat 
obscured by being in the wings of the uncoupled ab- 
sorption. In contrast, case (c) (A2 = —40 MHz, cor- 
responding to the minor dressed state) shows larger gain 
and exceptionally clear ringing. A fit to the decay of 
the ringing is shown in FIG. |^. As expected, this gives 
a decay rate of 5.5 MHz, which corresponds to half the 
phenomenological decay rate measured from the probe 
spectral linewidths. Finally, FIG. ||(d) is for a probe 
detuning of A2 = —29 MHz, again showing transient fea- 
tures between two regions of steady-state transparency. 



V. CONCLUSIONS 

We have carried out a theoretical and experimental 
study of transient phenomena associated with the switch- 
ing of the coupling field in a A type EIT system. These 
studies have shown novel transient gain features and aid 
the understanding of transient dressing. 

We have experimentally observed interesting transient 
gain features in both the turn-on and turn-off regimes. 
We have found it useful to interpret the transient dy- 
namics in terms of the filling and emptying of dressed 
states. The absolute peak gain values, especially the 
dressed state emptying, were found to be comparable 
in strength to the more usual, on-resonant turn-on case 
|p3| . Off resonant gain peaks have an important advan- 
tage over resonant gain peaks in that they can be tuned 
by appropriate choice of coupling strengths and detun- 
ings. In this way they can, for example, be removed from 
uncoupled absorptions or perhaps used to match transi- 
tions of different frequencies. The turn-off experiments 
also exhibited extremely clear transient signals and as 
such may be useful as a tool to explore the properties of 
atomic systems. 

We have obtained some new analytical results. In par- 
ticular we give the weak probe response after turning off 
the coupling field for arbitrary detunings. Analysis of 
the resonant turn-on using the Laplace method enabled 
us to identify a frequency oscillation not reported previ- 
ously. We analyzed this oscillation in terms of the 3D 
Vector Model and outlined a parameter regime where we 
would expect it to be observable. 

These results have applicability to the understanding 
of EIT for quantum information storage and may lead 
the way to methods which avoid the limitations of the 
condition of adiabaticity inherent in present schemes for 
optical quantum information storage j?], |f| . 
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APPENDIX A: FURTHER INVESTIGATIONS OF 
THE "FAST" OSCILLATIONS 

The oscillations in Im [/Ofc c (i)] at the angular frequency 
ill noted in Section [n] are unusual in that the Rabi fre- 
quency usually associated with the cycling of popula- 
tion is approximately f2i/2. The origins of this effect in 
the absence of decay are easily demonstrated using the 
3D Vector Model of 3-state systems introduced recently 
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The 3D Vector Model is based on the representation 
of the state of the atom \ip (t)) = v a \a) + Vb \b) — iv c |c) 
by the real, 3-dimensional vector lf(t) = (v a ,Vb,v c ) T , 
where here the states \a), \b) and |c) represent the bare 
energy eigenstates of the atom. The components of if 
are given with respect to three orthogonal axes which we 
shall call the a, b and c axes. These axes represent the 
three energy eigenstates of the atom. 

The action of the (resonant) probe and coupling fields 
on [?/;) is represented simply as a rotation about a Rabi 
vector fl = (—^fl 2 , ^f2i,0) T as follows: 



all? 



= n x 



(Al) 



The rotation of if about fl occurs at the frequency fl = 
| fl | = (fl'f + fi^) 1 ^ 2 /2 which is approximately fli/2 for 
f2 2 The Rabi vector lies in the a — b plane at an 

angle of 9 = tan^^f^/^i) ~ 0,2 /fli to the b axis. The 
value of pbc is the product of the amplitude of the state 
\b) and the complex conjugate of the amplitude of state 
|c) in \ip), which corresponds to the product ivbv c , and so 
Im(pf, c ) = VbV c . Absorption of the probe occurs when 
VbV c < and gain when VbV c > 0. 

We note from Eq. ( |l8| ) that, in general, Im (pbc) has fre- 
quency components at both the fli/2 and fli frequencies. 
However, there are two important cases where Im (pbc) 
oscillates at either fli or fli/2. 

The first case is when p® b — 1. Substituting this into 



Eq. (18) gives Im (pbc) oscillating at a frequency of fli/2. 
In this case the initial state vector lies along the ±b axis, 
that is, lf(0) — (0, ±1,0) T . In time the state vector if 
traces out an acute cone (of half angle 9) centered on 
± fl and revolves at frequency fl. The full solution of 
Eq. (Al) is given by 



(t)=± 



- cos 9 sin 9 [cos (fit) — 1} 
cos 2 9 + sin 2 9 cos (fit) 
— sin 9 sin (fit) 



(A2) 



Keeping terms up to first order in f^/fii gives 
lm[p bc (t)} w -(02/f2i)sin(Ot) ~ -(0 2 /Oi) sin(0 1 t/2) 
as found in Eq. (|i~8|). 



The second case is when p° aa = p bb = 1/2. In this case 
Eq. ( p"8| ) shows that Im (pb c ) oscillates only at the double 
frequency, i.e. fli. Equation (^) was derived for zero 
ground state coherence, i.e. p" b = 0. This means that in 
the gas sample, half the atoms are in state \a) and the 
remainder are in |6) . We model this by taking the average 
of the values of Im (pbc) for both collections of atoms. For 
the atoms beginning in the \a) state, the vector is initially 
aligned along the ±a axis, if (0) = (±1,0, 0) T , and traces 
out a cone of half angle ir/2±9. The solution of Eq. (Al) 
for this case is given by 



1 + [cos (fit) - 1] cos 2 9 
lf(t) = ±| cos 9 sin 9 [cos (fit) - 1] ). (A3) 
— cos 9 sin (fit) 

Keeping terms up to first order in f2 2 /^i gives 
lm[p bc (t)} w -(fi 2 AM [isin(2fii) -sin(fit)]. Tak- 
ing the average of this value and the value in 

Eq. (|A|) yields Im p bc (t) w -(fl 2 /fli)\ sin(2fii) « 
-(0 2 /r2i)|sin(fiit) as found in Eq. (|l|). 

In summary, we have traced the origin of the oscilla- 
tions at twice the Rabi frequency fl\ to the fact that for a 
pure state, pb c is the product of the amplitude of finding 
the atom in state \b) and the complex conjugate of the 
amplitude in state |c). In the absence of decay, both am- 
plitudes have terms that oscillate at the frequency fli/2. 
and so their product gives rise to a term oscillating at 
twice this frequency. So even though the populations os- 
cillate at the Rabi frequency fii/2, the absorption and 
emission of light has a term which occurs at twice this 
frequency. 
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